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ABSTRACT 

We co nstruct lightcones for the semi-analytic galaxy formation simulation of lGuo et al.l 
(|201ll) and make mock catalogues for comparison with deep high-redshift surveys. Pho- 
tometric properties are ca l culated with tw o different stellar population synthesis codes 
(|Bruzual fc Charlotll2003t lMarastonll2005l ) in order to study sensitivity to this aspect 
of the modelling. The catalogues are publicly available and include photometry for a 
large number of observed bands from 4000A to 6/xm, as well as rest-frame photom- 
etry and other intrinsic properties of the galaxies (e.g positions, peculiar velocities, 
stellar masses, si zes, morph o fogies , gas fractions, star formation rates, metallicities, 
halo properties). iGuo et all ()2011|) tuned their model to fit the low-redshift galaxy 
population but noted that at ^ 1 it overpredicts the abundance of galaxies below 
the "knee" of the stellar mass function. Here we extend the comparison to deep galaxy 
counts in the B, i, J, K and IRAC 3.6/^m, 4.5/xm and 5.8/im bands, to the redshift 
distributions of K and 5.8/im selected galaxies, and to the evolution of rest-frame lu- 
minosity functions in the B and K bands. The B, i and J counts are well reproduced, 
but at longer wavelengths the overabundant high-redshift galaxies produce excess faint 
counts. At bright magnitudes, counts in the IRAC bands are underpredicted, reflecting 
overly low stellar metallicities and the neglect of PAH emission. The predicted red- 
shift distributions for K and 5.8/im selected samples highlight the effect of emission 
from thermally pulsing AGB stars. The full treatment of Maraston (2005) predicts 
th ree times as many z ^ 2 galaxies in faint 5.8/im selected samples as the model 
of iBruzual fc CharlotI ()2003D . whereas the two models give similar predictions for K- 
band selected samples. Although luminosity functions are adequately reproduced out 
to z ^ 3 in rest-frame B, the same is true at rest- frame K only if TP- AGB emission is 
included, and then only at high luminosity. Fainter than L^, the two synthesis models 
agree but overpredict the number of galaxies, another reflection of the overabundance 
of ^ lO^^M© model galaxies at z ^ 1. 

Key words: methods; numerical - methods: statistical - galaxies: formation - galax- 
ies: evolution - stars: AGB 



1 INTRODUCTION 

Semi-analytic models of galaxy formation aim to predict 
the evolution of population properties such as the distri- 
butions of stellar mass, luminosity, star formation rate, 
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size, rotation velocity, morphology, gas content and metal- 
licity, as well as the scaling relations linking these prop- 
erties. They follow astrophysical processes affecting the 
baryonic components using a series of analytic, physically 
based model s which are embedded either in an analyt i c rep- 
resent ation (White & Frcnk 1991; Kauffm ann et al.l Il993l : 
ICole et al..,1994 : .Somerville fc Primack. 19991 ') or in a direct 
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numerical s imulation ijKauffmann et al ]| 19991 : ISpringeT et alj 
I2OOII . I2OO5I ) of the evolution of the underlying dark mat- 
ter distribution. Uncertain efficiencies and scalings of these 
astrophysical processes are represented by adjustable pa- 
rameters. These may be set a priori through a detailed 
calculation or simulation of specific processes, or they 



chosen data (e.f 


t. Croton et al.' 2006: Bower ct al.' 20061: 


Menci et al.ll2006 


:ICattaneo et al.. 2006: Dc Lucia & Blaizot 


2007 


: Monaco et al.l 20071: Somerville et al.ll2008l: IGuo et al. 


2011 


). The extremely broad range of relevant data and the 



considerable freedom in specifying appropriate recipes com- 
plicate the systematic comparison of semi-analytic mod- 
els with data. New and robust statistical tools have 
recently been developed to facilitate quantitative com- 



parisons ('Kampakoglou et al 
Henriques & Thomas 
20ld)- 



2010; 



20081: i Henrique s et al.l l2009l: 



Bower e t al. 2010l: ILu et al.1 



Such comparisons are sensitive to stellar population 
synthesis models which are required both to derive in- 
trinsic galaxy properties such as mass, age and star for- 
mation rate from observational data, and to calculate lu- 
minosities, colours and spectra for model galaxies. Erro- 
neous conversions between physical and observable prop- 
erties lead to incorrect conclusions about galaxy forma- 
tion physics, so it is important to check the implica- 
tions of ad opting differing stellar populations synthesis 



tions 01 ad optmg dittermg stellar populations synthesis 
models fe.g. iBuzzonil [19891 IWorthe^l 19941: IVazdekis et all 
19961: iFioc fc Rocca-Volrnera ngjll99?:^Leither er et al.lll999l: 
Bruzual fc CharlotI l2003l : iThomas ct al. 20031: iMarastonI 
20051 : IConrov et al.ll2009l ). 



Bruzual fc CharlotI l2003l : iThomas ct al. 20031: IMarastonI 
20051 : IConrov et al.ll2009l) . 

For example, the impact of including models for 
thermally-pulsating asymptotic giant branch (TP-AGB) 
stars has been studied in some detail in recent years. The 
contribution from these stars can significantly enhance the 
near- infrared emission of galaxies with Gyr-old populations 
(e.g. IMarastonI ITggi: iMarastoiil l2005l: Ivan der Wei et al.l 



2006I. iMarigQ fc Girardil l2007l : ICharlot fc Bruzual! [20071 : 
Conroy et al."200^. The data of IConrov et al.1 (|2009| ) and 



Zibetti e t al. (2009) show that inclusion of this additional 
emission can reduce the masses inferred from K-band 
light by as much as 0.6 d ex. The semi-analyti c mod els of 
Tonini et all (l2009l. l2010h . iFontanot fc Monacd (|2010l ) and 



Henriques et al.l ~ 2010^ "suggest that a substantial contri- 



bution from TP -AGB stars, as predicted by the model of 
IMarastonI (|2005h may explain the large number of extremely 
red objects found at 2 ^ 2. Other examples of how uncer- 
tainties in stellar evolution modelling affec t phys ical infer- 
ences from data are given bvlConrov et al.| |2010h . Here we 
compare predictions of the IGuo et aL ( 201 ll ) semi-analvtic 
model for two different population synth esis models, the one 
originally used by these authors (from iBruzual fc CharlotI 
[2003) and that of Maraston (20()1). 

Even neglecting uncertainties from stellar population 
modelling, mass-to-light ratios and other physical properties 
are often poorly constrained by available data. Estimates 
rely on fitting theoretical models to observed photometry 
and spectral energy distributions (SEDs) and approximately 
equivalent fits can often be obtained for broad ranges of 
assumed star formation history, chemical enrichment history 
and obscuration by dust. Additional uncertainties arise from 
possible variations in the initial mass function (IMF) with 



which stars form, and from possible spectral contributions 
from an active galactic nucleus (AGN). It has been argued 
in the past that a single optical colour (e..g. g — i) is, in 
practice, sufficient to deriv e light-to-mass ratios accurate to 
0.1 d ex for mos t galaxies (iBell et al.ll2003l : iGaUazzi fc Belli 
l2009h . However. lZibetti et al.l (|2009l ') showed that this is only 
true for relatively weak obscuration. For heavily obscured 
young populations, resolved photometry is needed to achieve 
an accuracy better than 0.2 dex and even that requires an 
additional near-infrared colour (e.g. one may combine g — i 
and i — H). 

Galaxy formation models directly predict star forma- 
tion and enrichment histories, so in the absence of obscura- 
tion a well defined SED can be predicted for each galaxy as 
a superposition of simple stellar populations (SSPs), each 
made up of coeval stars of a single metallicity. The "ob- 
served" photometry is then easily obtained by redshifting 
the SED and integrating over the appropriate photometric 
filter functions. In practice, however, the conversion to ob- 
servables is heavily influenced by dust and is sensitive to the 
detail s of its distributio n within a galaxy (e.g. lCranato et aP 
200C|; ICole et al.1 bood ) . This significantly hmits the preci- 
sion with which observables can be predicted from galaxy 
formation models. Current semi-analytic models often at- 
tempt to handle these uncertainties b y using observational 
data to constrain the dust model (e.g. Granato ot al. 200^ 
Cole et al.l I2OO0I : iKitzbichler fc Whitel ,2007: Guo fc Whit 



2009t ). 



Because of such difficulties it seems wise to compare 
theory and observation for a broad range of properties, at 
different redshifts, and at different "conversion levels". The 
latter is particularly crucial at high-redshift, where very lim- 
ited data are available and the relations between mass, light 
and star-formation rate are very uncertain. There is no pre- 
ferred "comparison frame" and conclusions are convincing 
only if a consistent picture emerges which matches srnoothly 
onto the lower redshift galaxy populations. IGuo et al.l(|201ll ) 
compare their model extensively to low-redshift galaxies but 
only present limited predictions at high redshift. Specifically, 
they compare to published estimates of the evolution of the 
stellar mass function of galaxies out to 2 ~ 4, finding signif- 
icant discrepancies for stellar masses below 5 X 10^° Mq. 

In this paper we extend this comparison considerably, 
analyzing the photometric properties of galaxies from high 
redshift to the present day, and comparing with observations 
at a variety of levels from number counts as a function of 
apparent magnitude, through redshift distributions of mag- 
nitude limited samples, to rest-frame luminosity functions as 
a function of redshift. In particular, we study predictions for 
the near-infrared bands for which data have recently become 
available from the Spitzer satellite. To facilitat e this work we 
build lightcones using the MoMaF software ijBlaizot et al.l 
[2005). We are then able to reproduce the observational selec- 
tion criteria of modern surveys such as the UKIRT In frared 
Deep Sky Survey (UKIDSS) llLawrence et al.l l2007h. th e 

t al.ll2 



20051 ). 



VIMOS-VLT Deep Survey (VVDS) l|Le Fevre et 
the Deep Evolutionary Exploratory Probe 2 Galaxy Redshift 
Survey (DEEP2) 



iry Jr'r 

(Davi s et al.| | 2003|) or the Cosmic Evolu- 



tion Survey (COSMOS) l|Scoville et allbOOTl). 

Si milar studies we r e perf orme d bv IKitzbichler fc Whitel 

(I2OO7I ), iGuo fc Whitl (|2009D and Ide la Torre et all (|201lD 

for earlier versions of the Munich semi-analytic model. The 
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results here are based on the model of IGuo et al] (|201ll ) 
which is implemented simultaneously on the Millennium 
and Millennium-II Simulations and was retuned to fit a 
broad range of "high precision" data on the low-redshift 
galaxy population, primarily from the Sloan Digital Sky 
Survey (SDSS). We also expand the photometric coverage 
from the ultraviolet to IRAC bands and test the depen- 
dence on stellar population synth esis modelling ov e r this 
wavelength range. In a recent paper, ISomerville et all (|201ll ') 
compared a different semi-analytic model with photometry 
extending to even longer wavelengths (the far infrared). This 
required modelling the re -emission of starligh t by heated 
dust, as also consider e d by Granato et al.l ||200(]|). I Cole et al] 
(|2000l ). iBaugh et all l|2005h , iLacev et all l|201ol ). Here we 
avoid this complication and compare only to directly ob- 
served starlight. 

This paper is organized as follows. In Section[2]we sum- 
marize the characteristics of the semi-analytic model we use 
and we describe how we construct lightcones for it. Section 
[3] then presents results for number counts and redshift dis- 
tributions as a function of apparent magnitude, and for rest- 
frame B- and if-band luminosity functions as a function of 
redshift. In Section [3] we present our conclusions. 



2 THE SEMI-ANALYTIC MODEL 

Modern semi-analytic models are built on merger trees 
from high-resolution dark matter simulations. These pro- 
vide a description of the evolution of the mass and num- 
ber density of dark matter halos and the subhalos within 
them, as well as of their spatial and kinematic distributions. 
The evolution of the baryonic components hosted by these 
(sub) halos is then followed using a set of simplified formu- 
lae describing each of the relevant astrophysical processes. 
The latest version of the Munich model (|Guo et al.i i2011] ) 
is implemented on two very large dark matter simula tions, 
the Millennium Simulation (MS; ISpri ngel et al.^ and 
the M illennium-II Simulation (MS-II; Bovlan-Kolc hin et al.l 
I2OO9I ). The MS follows the evolution of structure within a 
cube of side 500/i~'^Mpc (comoving) and its merger trees 
are complete for subhalos above a mass resolution limit 
of 1.7 X lO"/i'^M0. The MS-II follows a cube of side 
100/i~^Mpc but with 125 times better mass resolution (sub- 
halo masses greater than 1.4 x 10^ ft~^Mf7)). Both ado pt the 
same WM API-based cosmology jSpergel et aklbooj ) with 
parameters h — 0.73, fi,„ = 0.25, JIa = 0.75, n = 1 and 
as = 0.9. These are outside the region preferred by more 
recent analyses (in particular, erg appears too high) but this 
is of no consequence for the issues we study in this paper. 
For consistency, we will use this cosmology whenever it is 
necessary to derive the physical properties of galaxies from 
observed fluxes and redshifts. The distributions of physical 
properties converge in the two simulations for galaxies with 
10^'^Mq < M* < IO^ '^Mq. In this study we focus only on 
results from the MS, since its resolution limit is well below 
the stellar masses covered by the datasets with which we 
compare. 



2.1 The model of Guo et al. 2011 

For a full description of the se mi-analytic mod el used in 
this work we refer the reader to lGuo et all (|201ll ). Here we 
briefly describe changes from earlier versions of the Munich 
semi-analytic mod el that significantly a ffect our re sults. 



FoUowing iKitzbichler fc Whit j (120071') and 



iGuo fc Whiti ||2009D . the model of IGuo et al.l (|201lD 



includes a redshift-dependent model for internal extinction 
which assumes that the dust-to-gas ratio increases with 
metallicity but decreases with redshift. The effective optical 
depth is given by: 



10^^ atms cm 



r), (1) 



where (Nh) represents the mean column density of hydro- 
gen, {A\/Av)zq is the extinctio n curve for the solar metal- 
licity taken from iMathis et al.l ([l983) and s = 1.35 for 
A < 2000 A and s = 1.6 for A > 2000 A. 

When they implemented the lDe Lucia fc Blaizoj (|2007l ) 
version of the Munich model on the high-resolution 
Millennium-II Simulation. IGuo et al.l l|201l[ ) found it to over- 
produce dwarf galaxies. The authors therefore increased the 
efficiency of supernova feedback by introducing a direct de- 
pendence of the amount of gas reheated and ejected on the 
virial mass of the host halo. However, although the resulting 
model fits the stellar mass function of galaxies well at low 
redshift, it still produces more low-mass galaxies than are 
observed at z > 1. This deficiency is refiected in our results 
below. 

Finally, IGuo et all l|201lf ) introduced a more realistic 
treatment of satellite galaxy evolution and of mergers. The 
hot gas content of satellite galaxies is gradually stripped 
instead of being instantaneo usly removed at infall, a s sug- 
gested by the simulations of iMcCarthv et al.l (|2008l ). This 
allows satellites to continue forming stars for a longer pe- 
riod and reduces the excessively rapid reddening of these 
objects. In addition, satellites of satellites remain connected 
to their parent galaxies and can merge with them, rather 
than being automatically reassigned to the central galaxy of 
the group or cluster. The model also includes a treatment 
of the tidal disruption of satellite galaxies^ 



2.2 Lightcone Construction 

At high redshift, the observed fiuxes at a limited number of 
wavelengths are often the only data available for a galaxy, 
so that its redshift must be inferred through comparison of 
the observed colours to model templates. Even rest-frame 
magnitudes, colours and luminosities can then be subject 
to substantial uncertainties, and the conversion to intrinsic 
properties such as masses and star formation rates is prob- 
lematic. Results not only depend on the accuracy of the 
photometric redshift, but are also (almost) degenerate with 
respect to the star formation history, metallicity and dust 
content of the galaxy. These quantities are direct predictions 



^ See lHenriguea fc Thoma3 1I2OIC1I ) for an alternative extension of 
the Munich semi-analytic model modifying supernova feedback 
and including tidal disruption of satellites. 
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of a semi-analytic model, so that the conversion from intrin- 
sic to observed properties is, in principle, well defined, given 
a stellar population synthesis model, an assumed IMF and a 
specific model for intrinsic obscuration. It is thus often con- 
venient to consider conversion uncertainties as part of the 
model and to compare theoretical predictions directly with 
observables. To do this, we construct lightcones which allow 
the models to be "observed" in a way that mimics real sur- 
veys as directly as possible. We use two different population 
synthesis models to predict observables in order to assess the 
impact of the differing mass-to-light conversions they imply. 
We compute observed- and rest-frame fluxes from the ultra- 
violet to the near infrared so that our theoretical datasets 
resemble those of modern observational surveys not only in 
volume, but also in wavelength coverage. 

Our lightcones are bui lt using the Mock Map Facil- 
ity (MoMaF) developed bv iBlaizot et all (|2005h . We refer 
the reader to the original paper for a full description of the 
method. Here we briefly summarize the complications that 
arise when building lightcones from simulations of limited 
size and resolution. 

The Millennium Simulation has side of 500/i^^ Mpc 
(comoving). This is considerably smaller than, for example, 
the comoving distance to a galaxy observed at z ~ 2. Peri- 
odic replication of the simulation can lead to multiple ap- 
pearances of the same object within the lightcone, although 
typically at different redshifts and so with different prop- 
erties and at offset positions (due to large-scale motions). 
iBlaizot et al] (|2005h suggested applying a series of transfor- 
mations (rotations, translations and inversions) when tiling 
space with periodic replications. This does not, of course, 
prevent multiple appearances of a given object within the 
lightcone, but these duplicates are then viewed from differ- 
ent directions and no longer fall on a (nearly) regular lattice. 

Unfortunately, this technique also introduces disconti- 
nuities in large-scale structure at the boundaries between 
replications, affecting clustering statistics in a way which is 
at l east as difticult to model as that of the original periodic- 
ity. iKitzbichler fc Whitj (120071 ) showed that for lightcones 
of relatively small solid angle, the central line-of-sight can 
be chosen to pass through the lattice of periodic replica- 
tions in such a direction that multiple images of the same 
object are minimised or eliminated altogether. The latter is 
not possible if the comoving volume of the lightcone exceeds 
that of the simulation, but this technique can still be used to 
ensure that multiple appearances occur as far apart as pos- 
sible both on the sky and in redshift. We therefore use the 
method of IKitzbichler fc Whitd l|2007f ) in this paper. Space 
is flUed with periodic replications of the simulation, a posi- 
tion is chosen for the observer, and the central line-of-sight 
of the survey held is given a previously chosen orientation. 
Galaxies whose positions intercept the lightcone are selected 
and their comoving distance is converted i nto a r edshift. 

As explained in IKitzbichler fc White! (|2007l '). the time 
between stored snapshots for the Millennium Simulation 
varies between 100 and 380 Myr. This means that the in- 
trinsic properties of galaxies are not generally available at 
the time corresponding to their comoving distance. Rather 
they must be taken from the stored snapshot which is closest 
in distance to their light-cone position. Hence, galaxies with 
comoving distance {Di + Di+i)/2 > D > {Dt + _Di_i)/2 
are assigned the physical properties stored at Zi. The re- 



sulting discontinuity in galaxy population properties, at the 
boundaries between snapshots, could be reduced by inter- 
polating, but this works poorly for positions and velocities 
since the output separation is comparable to orbital times 
within groups and clusters. Moreover, it is not straightfor- 
ward for other galaxy properties either since these change 
discontinuously on timescales shorter than the output spac- 
ing, for example t hrough m ergers and starbursts. We thus 
follow iKitzbichler fc Whitd (|2007i ) and do not attempt any 
interpolation. 

The apparent luminosities and colours of galaxies de- 
pend strongly on their redshifts through the conversion 
between rest- and observed-frame photometric bands and 
through the inverse square dependen ce of apparent lumi- 
nosity on distance. As explained in IKitzbichler fc Whitd 
I2OO2) the flnal redshift of the galaxy in the lightcone is not 
available at the time observed-frame luminosities are com- 
puted in the semi-analytic model. However, there will be 
two extreme redshifts that bracket it. We compute apparent 
observed-frame luminosities (for flxed intrinsic properties) 
using these upper and lower limits, and once the galaxy 
is placed in the lightcone, we interpolate to obtain flnal 
observed-frame quantities. 

For this paper we construct lightcones for square areas 
of 1.4 X 1.4 deg^ out to high redshift with no faint magnitude 
cut. They are however limited by the mass resolution of the 
dark matter simulation (1.7 x 10^''/i~^Mq in halo mass) 
corresponding to stellar masses of ~ 10^-^Mq at z=0. While 
this does not matter for the questions we study in this paper, 
it should be borne in mind if the lightcones are used for other 
purposes^ 

2.3 Stellar Populations and Photometry 

2.3.1 Stellar Population Synthesis Models 

Semi-analytic models predict intrinsic properties of galaxies, 
such as stellar mass, star formation history, gas and dust 
content and metallicity. In order to convert these into ob- 
served spectral energy distributions (SED) or broad-band 
photometry, evolutionary population synthesis and dust 
models are required. The former predict the evolution of the 
light associated with a single short burst of star formation 
of given metallicity and with an assumed Initial Mass Func- 
tion (IMF), a so-called Simple Stellar Population (SSP). The 
intrinsic stellar emission from a model galaxy is then rep- 
resented as a superposition of SSPs weighted according to 
its star formation history. This emission must be processed 
through a dust model in order to predict the observable stel- 
lar emission. Uncertainties in the conversion between mass 
and light can jeopardise any comparison between theory 
and observations. There are still signiflcant differences be- 
tween published evolutionary population synthesis models 
and these should be considered as part of the systematic un- 
certainties when comparing semi-analytic model predictions 
to data. Throughout this paper we present results from two 
distinct stellar population synthesis codes : one that has b een 
traditionally used in the Munich model (jBruzual fc CharlotI 

^ With the acceptance of this paper the lightcones are publicly 
released at http:/ /www.mpa-garching.mpg.de/millenniumj 
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Figure 1. The response functions of tiie filters for which fluxes are 
computed for the lightcones produced in this paper. These extend 
from the far ultra-violet to the IRAC bands and include: in the 
top panel, GALEX FUV and NUV, Johnson UBVR^IcJHKsK 
and IRAC bands; in the second panel, SDSS ugriz; in the third 
panel, HST WFC3 UV and IR and ACS WFC; and in the bottom 
panel, HST WFPC2, VIMOS U and NICMOS. 



l2003h and the iMarastmJ fe005h mode l impl emented in the 
semi-analytic code by Hcn rinues et al.l l|2010l ). In both cases 
we adopt the same lChabrieij (|2003l ) IMF and a similar metal- 



licity grid. We hope that the differences we find will give 
some indication of the impact of mass-to-light conversion 
uncertainties on galaxy formation modelling. 



2.3.2 Photometry 

In order to increase the predictive power of the model and 
allow it to be tested against a wider range of observations, 
we also expand the number of photometric bands for which 
fluxes are computed, covering all wavelengths dominated by 
direct emission from stars, from the UV to the near-infrared 
IRAC bands. In Fig. [1] we plot the relevant filter transmis- 
sion curves. In the top panel we show the GALEX FUV and 
NUV, the Johnson U, B, V, Rc,Ic, Z, Y, J, H,K^,K and the 
IRAC 3.6/im, 4.5/xm, 5.8^m and S.Ofim bands; in the second 
panel, the SDSS u,g,r,i, z bands; in the third panel, bands 
from HST instruments, three UV bands from the WFC3- 
UVIS (0.225/im, 0.275/im, 0.336/im), seven optical bands 
from the ACS- WFC (0.435/im, 0.475^im, 0.606^m, 0.625/im, 
0.775/-im, 0.814/-im, 0.850/im) and three near-infrared bands 
from the WFC3-IR (1.05^im, 1.25Atm, 1.60^m); and in the 
bottom panel, the VIMOS U band, the 2 NICMOS near- 
infrared bands (1.1/im and l.Gfim) and two WFPC2 bands 
(0.30/xm and 0.45^m). While we will not give results for all 
these bands in this paper, we will include the relevant ap- 



parent magnitudes in our light-cone catalogues in order to 
enhance their utility to others. 

All magnitudes are in the AB system. In order to be as 
close as possible to observations, we use the Ks-ha,nd when 
presenting results for number counts and redshift distribu- 
tions and the .ft'-band when discussing the evolution of the 
rest-frame luminosity function. 

The lightcones constructed and made public in this 
work provide a useful tool to test observational derivations 
of intrinsic galaxy properties. The wide wavelength cover- 
age of observed- and rest-frame photometry, together with 
the two stellar population synthesis models considered, can 
be used to check derivations of rest-frame magnitudes from 
observed photometry, as well as the reliability of properties 
obtained from SED fitting, such as stellar masses, ages and 
star-formation histories. 



3 RESULTS 

In this section we compare predictions of our models to ob- 
servational data. We start with number counts as a function 
of apparent magnitude in a wide range of photometric bands 
(from the optical blue to the IRAC bands) and move on to 
redshift distributions for K and IRAC 5.8/-im selected galax- 
ies. Finally we investigate the evolution of the rest-frame op- 
tical and near-infrared luminosity functions which, although 
further from the directly observed quantities, allow a better 

und erstanding of galaxy e v olution . 

_ Kitzbichler fc White! (|2007l ') and Ide la Torre et all 



presented similar tests for an earlier version of the 
Munich semi-analytic model. Here, we take advantage of 
recent advances in the available observations and extend 
these comparisons to higher redshift and to a wider range of 
wavelengths. We also test the impact of population synthe- 
s is mo dels by comparing results for the iBruzual fc CharlotI 
(i2003j ) an d the Maraston (.2005.) m o dels. This follows up 
work by iTonini et al. "liogd 120101) : iFontanot fc Monaco! 
(|2010! ) and !Henriaues et al.! (!2010! '). who showed that the 
inclusion of near-infrared emission from TP-AGB stars 
increases the predicted number of massive and extremely 
red objects at z ~ 2, as seems to be required by observation. 
Our comparison is based on a large number of lightcone 
realizations with areas and selection effects matching the 
relevant observational surveys. 



3.1 Number Counts 

Galaxy counts in a given observed band can be difficult to in- 
terpret. At each apparent magnitude they consist of galaxies 
at a wide range of redshifts and thus with correspondingly 
wide ranges of absolute magnitude and of emitted wave- 
length. Nevertheless, such counts provide an important test 
of models because they are directly observed and so are in- 
dependent of uncertainties in redshift, k-correction, obscu- 
ration correction, etc. 

Fig. [21 shows galaxy number counts for the B and 
i, J, Ks, IRAC 3.6/im, 4.5/-im and 5.8/^m bands (from 
top left to bottom right). Solid red and dashed blue 
lines are model pr edictions for the iMarastonI (|2005! ) and 
iBruzual fc CharlotI (2003) stellar population synthesis mod- 
els respectively. Filled regions show the la field-to-field scat- 
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Figure 2. Galaxy number counts as a function of apparent magnitude. From top left to the bot tom right, t he pa nels s how number counts 
in the B and j, J, Ka, IRAC 3.6/jm, 4.5/jm and 5.8/xm bands. Theoretical predictions for the lMarastonI | |2005| ) and'Bruzu al fc Charloj 
1I2OO3I ) population synthesis models are shown as solid red and dashed blue lines respectively. The filled regions represent the Itr field-to- 
field scatter for surveys of area 2 deg^, except for the IRAC bands a t faint ma gnitudes (M>18.5) where 100 arcmin^ fields are assumed. 
The B-band numb er counts are compared with data from SDSS l lYasuda et al. 2001), VVDS (McCrackcn ct al. 2003 ) and HDF-N 
jCapak et al .120041) : i-band counts are also compared with data from SDSS a nd VVDS and with th e COSMOS sample jCapak et alj200^: 
for th e J- and Ka- b ands we s how observations from the CDF and HDF-S iSaracco et al.ll200ll ). DEEP2 and Palomar JCoii ^^^^^^^^' 



200S|) and MOIRCS jKeenan e t al. 2010) with UKIDSS-U DF data jCirasuolo et al.ll2010l) plotted for the iiT-band; Spitzer l lFazio et al. 



20041) , FIREWORKS llWuvts et al...200a) and NEWFIRM jWhitaker et al.ll201 j) data are shown for the IRACa.s^m, IRAC4.5^m and 



IRAC5 .8A[m bands. 



ter expected among surveys of area 2 deg , except that 100 
arcmin^ fields are assumed for the IRAC bands at faint mag- 
nitudes (M>18.5). 

The optical B- and i-band number counts are com- 
pared with data fro m the Sloan Digital Sky Survey (SSDS; 



lYasuda et al. 2001 ) and the Vimos Very Deep Survey 
fVVDS: lMcCracken et al.ll2003l ') at brighter magnitudes. At 
fainter magnitudes we use the Hubble Deep Field - North 
(HDF-N; Cap ak ,ct_al. 2004) for the B-band and the COS- 
MOS sample (|Capak et~alll2007l ') for the i-band. Both pop- 
ulation synthesis models match the data for bright galaxies 
(seen at low redsh i ft in the rest-frame optical), while the 
iBruzual fc Charlotj l| 20031 ) model predicts more galaxies at 
faint apparent magnitudes, in better agreement with obser- 
vations. Similar trends were fou nd for previous versions of 
the Munich semi-analy tic model (jKitzbichler fc Whit3l2007l : 
Ide la Torre et al.ll201ll ). It is difhcult to draw firm conclu- 
sions from this result, however, since these faint counts cor- 



respond to rest-frame ultra-violet emission from galaxies at 
high redshift where uncertainties affect not only the stel- 
lar population modelling, but also the simplistic treatments 
of starbursts and of dust o bscuration in the semi-analytic 
model. Conrov et al.l (|2010l ) showed, for example, that in- 
creasing the number of blue stragglers or blue horizontal 
branch stars increases the predicted ultra-violet emission 
from passive galaxies. 

We compare the J- and ATs-band counts to observations 
from the Chandra Deep Fie ld and the Hubble D eep Field 
South (CDF and HDF-S; ISaracco et all I2OOI]) fro m the 



DEEP2 and Palomar surve ys llConselice et al. 20081 ) . and 
from the MOIRCS sample (iKeenan et all I2OI0I ) . In addi- 
tion, we show K-band counts b ased on the UKIDSS -Ultra 
Deep Field data (UKIDSS-UDF: ICirasuolo et al.ll2010l ). The 
two stellar population synthesis models give similar predic- 
tions for the J-band number counts which agree with the 
data. Both predict too many faint objects in the Ka and K 
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bands. As shown bv lGuo et aLll|201ll ). and as this paper will 
clarify, this is because the semi-analytic model overpredicts 
the abundance of low-mass galaxies at high redshift. The two 
stellar population synthesis models predict similar counts at 
both bright (low redshift, near infrared emission) and faint 
(high redshift, red optical emission) apparent magnitudes, 
but they disagree at intermediate apparent magnitudes. As 
we will see in more detail below, the difference is a con- 
sequence of TB-AGB emission from stars with ages of one 
or two Gyr w hich is fully included in th e iMarastonI (|2005l ) 
but not in the iBruzual fc Charloj l|2003l ) stellar population 
model. 

Predicted number counts for the IRAC S. G^tm, i.bfim 
and 5.8pm bands are plotted against Sp itzer (|Fazio et al.l 
|2004). FIREWORKS l|Wuvts et all |2008| ) and NEWFIRM 
l|Whitaker et al.ll201ll ) observations. Both the models and 
the observations show a pronounced change in slope at an 
apparent magnitude near 20, but the break is stronger in 
the observations than in the model and occurs at slightly 
brighter apparent magnitudes. As a result, the models 
under-predict the number of bright objects (low redshift, 
emission longwards of the rest-frame K-band) and over- 
predict the number of faint objects (high redshift, emission 
in the rest-frame JHK region). The latter underprediction 
is even more pronounced here than in the K-band and again 
is likely due to the overabundance of lower mass galaxies at 
a ^ 1 in the model. The deficit of bright galaxies is visi- 
ble also in the z = Q rest-f r ame J f-band luminosity function 
(Fig. [5]). Since [Guo et al.l (I2OIII ) tuned their semi-analytic 
model to match the observed low-redshift stellar mass func- 
tion, this deficit implies overly large mass-to-near-infrared- 
light ratios which might b e explained by overly smal l 
stellar metallicities. Indeed, iHenrigues fc Thomas! (|2010| ) 
showed that the most ma ssive low-redshift galaxies in the 
Lucia fc Blaizod (|2007l ) version of the model have stellar 
metallicities which are too low by about a factor of two (the 
dashed red lines in their Figs 4 and 10) . An increase in metal- 
licity could remove the discrepancy by reducing the mass- 
to-near-infrared-light ratios in model|j Another important 
factor, particularly for the IRAC 5.8/im band, is the possible 
contamination by emission from hot dust, spec ifically, emis- 
sion i n the 3.3, 6.2 and 7.7nm PAH features JPraine et al] 
I2OO7I : iDraine fc Lill2007l : Ida Cunha et"aLll2008l ). Such emis- 
sion is not included in the models but may well be significant 
in the real low-redshift galaxies. 

3.2 Redshift distributions for K and IRACs.g^m 
selected samples 

In Fig. [3] we compare predictions from our models to the 
observed photometric redshift distributions of galaxy sam- 
ples selected above Kg and IRAC 5.8/xm apparent mag- 
nitude limits. As for the number counts, the solid red 
and dashed blu e lin es represent predict i ons b ased on the 
iMarastonI (120051 ) and lBruzual fc CharlotI (|2003l ') stellar pop- 
ulation synthesis models. Filled regions show the expected 

^ The luminosity correction depends on the actual deficit in 
metallicity, which in turn depends strongly on which stellar pop- 
ulation model is used to derive masses and metallicities for the 
observed galaxies. 



1(7 scatter among fields with an area of 100 arcmin^. In the 
left panel the number of galaxies per unit area and redshift 
is plotted for samples with Kg < 21.8 and Ks < 23.3, while 
the right panel gives similar results but for samples with 
IRAC 5.8/im apparent magnitude brighter than 21.8 (and 
Ks < 23.0). Our theoretical predictions are compared with 
data from two publ ic catalogues: the F IREWORKS data for 
the GOODS-CDF (|Wuvts et al.ll2008l l and the NEWFIRM 
Mediu m Band Survey data for the COSMOS and AEGIS 
fields jWhitaker et al.ll201ll ). The wide photometric cover- 
age of these two datasets results in robust and relatively 
precise photometric redshift measurements. Our theoretical 
samples are selected using photometric criteria very similar 
to those defining the observed samples, although we plot the 
distribution of their true redshifts rather than attempting to 
reproduce the observational redshift estimation procedure. 

Selection by observed-frame Ks-hand magnitude picks 
galaxies on the basis of their rest-frame -ft'-band emission 
at low redshift, their rest-frame J band emission at z ^ 1, 
and their rest-frame optical emission at z > 1.5. For both 
magnitude limits, the samples are dominated by intrinsi- 
cally faint objects at low-redshift but by intrinsically bright 
galaxies beyond z ~ 1 (for Ks < 21.8) or 2; ~ 1.5 (for 
Ks < 23.3)). For both apparent magnitude limits the pre- 
dictions for the two population synthesis models agree, and 
for the brighter limit they are consistent with the observa- 
tions out to a redshift of almost two. Both underpredict the 
counts at higher redshift, with th e effect being slightly larger 
for the lBruzual fc CharlotI (|2003l ) model. This may refiect an 
underabundance of intrinsically bright objects (in the opti- 
cal) in the model at these redshifts, but could also be due to 
magnitude and photometric redshift errors in the data which 
primarily affect the tails of the distribution. For the fainter 
apparent magnitude limit the model clearly overestimates 
the number of objects over the redshift range 1 < z < 2.5. 
These galaxies typically have stellar masses of order a few 
10^" M0 and this discrepancy refiects the overabundance of 
objects of t his m ass and redshift unity that was flagged by 
iGuo et al.l (|201ll '). The same problem was idcntifled in ear- 
lier versions of the model by [Kitzbich lcr fc White (200'i) 
and Ide la Torre et all (|201ll ). At higher and lower redshifts 
the abundances agree quite well in model and data, reflect- 
ing the fact that the semi-analytic model was tuned to flt 
galaxy abundances at low redshift, and predicts an abun- 
dance of high-mass galaxies which flts observed estimates 
quite well at high redshift. 

For the IRAC 5.8/^m selected samples shown in the 
right panel of Fig. [31 there is a signiflcant difference be- 
tween the pred ictions of the two p o pulat ion synthesis mod- 
els. While the iBruzual fc CharlotI l|2003l ') model predicts a 
dist ribution with sim ilar shape to those in the left panel, 
the IMarastonI (|2005l ) model makes a concordant prediction 
only at z < 1.2. Beyond this point there is a "bump" and at 
higher redshift it predicts roughly 3 times as many galaxies 
as the IBruzual fc CharlotI l|2003l ) model. A corresponding 
bump is not present in the observational data w hich are 
better described by the IBruzual fc Ch ariot' (2003") model, 
at least out to z ^ 2.5. The bump in the Maraston (200^) 
model is caused by strong rest-frame JHK emission from 
TP-AGB stars associated with intermediate age stellar pop- 
ulations. While this emission brings the predicted numbers 
of galaxies into rough agreement with the data at the high- 
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Figure 3. The redshift distribution of galaxies selected above observed-frame i^s-band (left panel) and IRAC 5.8/im-band (right panel) 
apparent magni tude limits. The solid red and dashed blue lines represent the mean predictions of our semi-analytic model for the 
iMarastonl l|2005l ) and lBruzual fc CharlotI l|2003l ) stellar population synthesis models. The filled regions show the la scatter among fields 
with area 100 arcmin'^. The left panel shows the total number of galaxies per unit area and redshift for samples with Ks < 21.8 and 
Ks < 23.3. Similar curves are shown in the right panel but for samples with IRACs.gfim <21.8 and Kg <23.0. Theoretical distributions 
are compared with data from FIREWORKS (,Wuvts ct al.,.2008 ^ and from the NEWFIRM Medium Band Survey ^Whitaker et al..,201ll ). 



est redshifts, it results in an overabundance at 2 ~ 2 where 
the observational datasets appear most robust. Since this 
effect is also present for data in other wavebands, for which 
TP-AGB emission is not an issue, it suggests that it results 
from the semi-analytic model overpredicting the abundance 
of the relevant moderate mass galaxies by a substantial fac- 
tor at this redshif t. In view of this, the agreement achieved 
by the |Bruzual fc CharlotI ([2003) model is probably coinci- 
dental, resulting from the overestimated abundance of mod- 
erate mass galaxies being compensated by an overestimate 
of their rest-frame JHK mass-to-light ratios. 

3.3 The rest-frame B-Band Luminosity Function 

Rest-frame luminosity functions and colour distributions 
as a function of redshift provide direct estimates of the 
abundance evolution of various galaxy types (e.g. star- 
forming/passive, high/low mass). However, they require ac- 
curate redshifts and appropriate photometry if they are to 
be determined reliably from observed-frame fluxes. The wide 
wavelength coverage of modern surveys produces robust 
photometric redshifts, and, in addition, allows rest-frame 
optical and near-infrared magnitudes to be determined by 
interpolation over the full range < 2 < 4, rather than 
req uiring an un c ertain extrapolation based on an SED fit. 

IGuo et al] (|201ll ) showed that their semi-analytic 
model reproduces observed z ~ 0.1 luminosity func- 
tions in the SDSS g,r,i and z bands. At higher redshift 
they implemented t h e red shift-dependent dust model of 
iKitzbichler fc White! (|2007l ). This reproduces the observed 



abundance of colour-selected galaxies at z ^ 2 and 2 ~ 3 
for the pr evious version of the semi-analytic model (see 
iGuo fc Wh ite 2009). In Fig. H 

we show the evolution of 
the B-band luminosity function from 2 = to 2 = 3 
for our current semi-analytic model. Solid red a nd da shed 
blue lines, represen t versi ons with the IMarastonl |200i) and 
iBruzual fc CharlotI (|2003h stellar population models respec- 
tively. Filled regions give the expected la field-to-field scat- 
ter for surveys of area 1.4 deg'^, except that smaller fields 
(with area 150 arcmin'^) were assumed for the 2.5 < 2 < 3.5 
panel, for the 1.3 < z < 2.0 panel fainter than —21.0 and 
for the 0.8 < z < 1.2 panel fainter than —19.0. 

At 2 ~ the model 6j-band luminosity functiorQis com- 
p ared with the 6 Deg ree Field Galaxy Redshift Survey result 
of iJones et al.l l|2006l . 6DFGRS), repeated for reference as a 
black dotted line in all panels. For the 2 ~ panel we use the 
final snapshot of the simulation rather than the lightcone, 
finding excellent agreement with the 6DFRS result, just as 
was the case for the corresponding SD SS luminosity func- 
tion (in the g band) in lCuo et all (|201lh . For the other 2^1 
panels we compare with data fro m the relatively wide V VDS 
(|llbert et all |2005|). DEEP2 jWillmer et all I2OO6I ) and 
zCOSMQS (|Zucca et aLll2009l ) surveys. At higher redshift 
only data for smaller fields ar e available. We show results 
fromPoli et al." ("2003!, HDF-S). [Gia,llongo et all l|2005l . HDF- 
N)jMarchcsini ct al. (20o3, GOODS-MUSYC-FFIRES) and 
ISalinAeni et all (20081 . GOODS-MUSYC). 



" We assume bj = B - 0.267(5 - V), l lNorbere et alj|2002t) . 
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Figure 4. Evolution of the r e st-fra me B-band luminosity function from 2 = 3 to 2 = 0. Theoretical predictions for the iMarastonl 
1I2OO5I ') and lBruzual fc Chariot! i200^ stellar population models are shown as solid red and dashed blue lines respectively. Filled regions 
represent the la field-to- field scatter expected for surveys of area 1.4 deg-^, except in the highest redshift panel, where 150 arcmin^ fields 
are assumed. Fields of this size are also assumed for the intrinsically fainter galaxies in the 0.8<z<1.2 and 1.3<z<2.0 panels (for galaxies 
with > —19.0 and > —21.0 respectively). At 2 = t he model bj-band luminosity function {bj = B — 0.267(i? — l/);^Norbcrg ct ahl 
I2OO2I) is compared with observations from the 6DFGRS lljones et aI.ll2006^ . repea ted at all redshifts a s a dotted bl ack line. At highe r 
redshift s we show obser-y ational est imates from VVDS llllbert et al.ll2005^ . DEEP 2 llWillmer et al.ll2006h. zCOSM OS (tZucca et al.ll2009l ). 
HDF-S llPoh et al.ll2003f) , HDF-N jCiallongo et al.ll2005l '), GOODS-MUSYC plus FIRES jMarchesini et ai]|2007^ and GOODS-MUSYC 
dSalimbeni et ahlbOO^T ^ 



The model reproduces the evolution of the rest-frame 
B-hand luminosity function reasonably well out to 2: = 3. 
It overpredicts the abundance of faint objects at 2: ~ 2 and 
underpredicts the abundance of bright objects at 2 ~ 3, al- 
though it may still be compatible with the data given the rel- 
atively large error bars quoted by the observers and the sub- 
stantial scatter between the observational determinations. 
The two stellar population synthesis models give very sim- 
ilar results in this band. We note that the predicted fluxes 
are strongly affected by dust, and so are dependent on the 
adopted dust model. Further testing of the simplistic and 
relativ ely poorly motivated model of iKitzbichler fc White! 
I2OO7I ) is clearly needed. 



3.4 The if-Band Luminosity Function 

The K-hand luminosity function has long been thought of 
as a proxy for the stellar mass function. Recent results 
have shown, however, that this assumption, while moder- 



ately accurate at low redshift, can break down badly at 
early times. Notably, the fact that the characteristic lumi- 
nosity L* incr eases with increasing redshift just as for the 
optical bands jCirasuolo et al.l!2010l ) is inconsistent with a 
time-independent K-hand mass-to-light ratio, which would 
imply the "evaporation" of material from the most massive 
galaxies. This luminosity function behaviour is easily under- 
stood in the context of recent stellar population synthesis 
models. In particular, a significant amount of if-band emis- 
sion comes not from the old populations which dominate 
the stellar mass, but rather from intermedi ate age stars (~1 
Gyr) passing through th e TB-AGB phase (!Marastonl!2005! : 
!Charlot fc Bruzuall!2007! ). At high redshifts these relatively 
young populations can dominate the rest-frame luminosity 
in the Ti'-band and they are only later replaced by predom- 
inantly old populations (|Henriaues et al. 2010). 

In Fig. [S] we plot the evolution of the Ji'-band luminos- 
ity function out to z = 3. As in previous figures, the solid 
red and dashed blue lines represent predictions based on 
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Figure 5 . Evolu tion of the rest-frame K - ha.nd luminosity function from z = 3 to z = 0. Predictions of our semi-analytic model for the 
iMarastonI | |2005|) and lBruzual fc Charlol ||2003| ') population synthesis models are shown by solid red and dashed blue lines respectively. 
Filled regions represent the expected l a field-to-field sc atter for surveys of area .7 deg;-^. At z "-^ t he model X^-band luminosity function 
is compared with data from 2MASS jBell et al.ll2003l'l and 6DFGRS+2MASS l ljones et aLll2006l). We repeat the l atter at all redshifts 
as a black dotted line . At higher redshi fts we show observat ional estimates b ased on MUNICS llDrorv et al.ll2003t) . the UKIDSS-UDF 
dCirasuolo et al.ll20l6l ). the K20 Survey fPo zzetti et al.ll2003l) and the HDF-S jSaracco et al .' 2006). Note that in all surveys other than 
the UKIDSS-UDF, the rest-frame K-luminosities are not directly measured but are rather estimated by extrapolating from the observed 
frame iC-band fluxes using an SED model. 



the IMarastonI (120051 ) and iBruzual fc Charloj l|2003h stellar 
population models, and filled regions outline the expected 
la field-to-field scatter among surveys of area 0.7 deg^. At 
2 ~ the model Ka-hanA luminosity function is comp ared 
with observational data from 2MASS (|BeU et allbood ) and 
6DFGRS-H2MASS (| Jones et al.ll2006l ). As a reference, we 
repeat the latter at all redshifts as a black dotted line. 
For the 2 ~ panel we use the final snapshot of the 
simulation rather than the lightcone to obtain the theo- 
retical prediction. At higher redshifts, the model is com - 
pare d with data from the MUNI CS jProrv et al.' 200^), 
K20 (iPozzetti et"al. 2003). HDF-S llSaracco et al. 2006) and 
UKIDSS-UDF lCirasuolo et all (|2010h . Only the last of these 
surveys estimates the rest-frame /C-band flux directly by 
interpolating between observed-frame magnitudes at corre- 
sponding wavelengths (from Spitzer/IRAC). The other sur- 
veys extrapolate the observed-frame K flux to longer wave- 
length using an uncertain SED fit and thus may be subject 
to substantial systematic errors. 

Our two population synthesis models give very similar 



predictions for the rest-frame K-hand luminosity function 
out to z ~ 0.5. At higher redshifts their shapes and their 
normalisations remain si milar but their ch aracteristic lumi- 
nosities diverge with the IMarastonI l|2005l ) prediction being 
brighter by about 0.4, 0.5, 0.6 and 0.75 magnitudes at red- 
shifts of 1.0, 1.5, 2.0 and 3.0 respectively. This refiects the 
increasing contribution from TP- AGE stars as the mean age 
of the galaxies gets younger. The predictions of both models 
are stron gly at variance w ith observation at these redshifts. 
While the lMarastonI l|2005l ) model agrees with the high-mass 
tail of the observed luminosity functions at all redshifts, it 
seriously overpredicts the abundance of less massive galax- 
ies aX z = 1 and earlier. For the lowest luminosity bin of 
the UKIDSS-UDF dataset the overprediction is by factors 
of 2, 4, 6 and al most 8 at redshifts o f 1.0, 1.5, 2.0 and 3.0 
respectively. The ISruzual fc CharlotI (|2003l ) model fails to 
reproduce the rest-frame K luminosities of the most mas- 
sive systems (by about 0.7 mag by z ~ 3) but nevertheless 
overpredicts the abundance of less luminous systems almost 
as badly as the iMarastoui ((2005i ) model. This is the sub- 
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stantial problem already pointed out bv lGuo et alj (|201l[ ): 
their galaxy formation assumptions produce moderate mass 
galaxies (M* ~ 10^° M©) too early to be compatible with 
current data on populations at z ^ 1. 

As noted above when discussing Fig. (2] the model 
also slightly underestimates the K-hand luminosities of 
massive low-redshift galaxies. This more subtle problem 
is due to these massive galaxies being too blue (s ee the 
colour distributions in Fig. 12 of IGuo et al.l l201lh . This 
is likely caused by aii unde rabundance of heavy elements 
l|Henriques fc Thomaslboid '). 

Similar results for the evolution of the rest-frame 
K-h&nd luminosity function have been obtained for 
an earlier version of th e Munich semi-analytic model 
l|De Lucia fc BlaizotI |2007| as w e ll as for the semi-an alyt ic 



models of 



Menci et alT 120061, [Monaco et"al] |2007| and 



iFontanot et al.l l2009l ). In a recent paper ISomerville et al.l 

compared another independent semi-analytic model 
to observational data on the evolution of the rest-frame 
1500A, B and K band luminosity functions. The authors 
were able to get a reasonable match to the bright tail with- 
out including TP-AGB emission, but they do obtain a much 
weaker evolution of the characteristic Lj, t han o bserved, 
in concordance with our iBruzual fc Chariot! l|2003l ) results. 
They also overpredict the abundance of lower luminosity 
galaxies by very similar factors to those that we find here. 
It seems that whatever is causing the overly early formation 
of lower mass galaxies is common to all recent semi-analytic 
models. 



4 CONCLUSIONS 

We have constructed lightcones fro m the latest ve rsion 
of the Munich semi-analytic model ||Guo et al.l 1201 ll ) and 
used them to compare the model with the high-redshift 
galaxy population as revealed by recent deep surveys at 
optical and near-infrared wavelengths. We have combined 
the mode l with two different stellar population sy nthesis 
packages (jBruzual fc Chariot! l2003l : iMarastonI I2OO5I ) in or- 
der to understand how differences in the photometric mod- 
elling are reflected in inferences about galaxy evolution. We 
use multiple independent lightcones to characterize cosmic 
variance uncertainties in currently available datasets. Our 
mock catalogues are made publicly available with publica- 
tion of this paper and provide observer-frame photometry in 
40 commonly used photometric bands, in addition to rest- 
frame photometry and a variety of physical properties of 
the galaxies (positions, peculiar velocities, stellar masses, 
halo masses, sizes, morphologies, gas fractions, star forma- 
tion rates, metallicities, halo properties). 

We now summarise the principal conclusions from our 
comparison of models and data. 

• For both stellar populations the model matches the 
observed- frame B, i and J number counts but overpredicts 
the counts at faint magnitudes {niAB > 20) in the Ks, 
IRAC 3.6/im, 4.5fim and 5.8fim bands. This reflects the 
overproduction of moderate mass galaxi es (stellar masses 
~ lO^^Mo) at 2 > 1 already noted bv lGuo et all l|201ir ). 
The matching of the faint optical counts is fortuitous - this 
overproduction is masked by overly large mass-to-light ra- 



tios in the rest-frame near-UV, perhaps due to problems 
with the dust modelling. 

• At bright magnitudes (ttias < 20) the model under- 
predicts the counts in the three IRAC bands. This is due 
to an underestimation of the near-infrared luminosities of 
low-redshift massive galaxies caused in part by the fact that 
such galaxies are insufficiently metal-rich in the model, and 
in part by the model's neglect of PAH emission from hot 
dust (which is particularly significant at 5.8/im). 

• At magnitudes where the model Ks-hand counts agree 
with observations the redshift distribution of A's-selected 
samples is also reproduced. At fainter magnitudes where the 
counts are overpredicted, the excess galaxies occur primar- 
ily at 1 < z < 2.5, again refiecting the overproduction of 
Mi, ~ 10^" M0 gala:xies at these epochs. The two population 
synthesis models give similar results in both regimes. 

• The two population synthesis models predict differ- 
ent redshift distributions for galaxies selected to ttias ~ 
22 in the IRAC 5.8^m band. Emission from TP-AGB 
stars enhances t he number of galaxies at z > 1.5 in the 
'Maraston (200^ mode l by about a factor of 3 relative to 
the Bruzual fc CharlotI (,2003) model, causing it to over- 
predict the observ e d abu ndance at z ~ 2. Although the 
iBruzual fc CharlotI (|2003h model agrees with the observed 
redshift distribution for < < 3, this is a result of the 
overabundance of moderate mass galaxies being cancelled 
by an overestimate of their near-infrared mass-to-light ra- 
tios. 

• The two population synthesis give similar results for 
the evolution of the rest-frame B-band luminosity function, 
agreeing well with observation out to z — 1.2. At higher red- 
shift the agreement is less convincing. The overabundance 
of lower mass model galaxies starts to become evident, and 
there is some indication that the models underpredict the 
abundance of the most luminous objects. Cosmic variance 
and other uncertainties in the currently available data, to- 
gether with dust modelling uncertainties in the model, pre- 
clude an y strong conclusi ons. 

• The IMarastonI (12005) population model reproduces the 
bright tail of the rest-frame J^-b and luminosity function al l 
the way out to z ~ 3, whereas the lBruzual fc Chariot! (|2003l ) 
model underpredicts the near-infrared luminosities of these 
massive galaxies by an amount which increases from about 
0.3 mag at 2 ~ 1 to 0.7 mag at z ~ 3. The overproduction of 
Af* ~ 10^" Af© galaxies at these times causes both models to 
substantially overpredict galaxy abundances below the knee 
of the luminosity function. 

In the literature it has often been suggested that 
semi-analytic models fail to reproduce the rest-frame K- 
band galaxy luminosities of the brightest high-redshift 
galaxies (at z ^ 2 — 3), and the failure is usu- 
ally attribute d to insufficie n tly r a pid mass growth at 
early times (iPozzetti et al.l [200^: ICimatti et all 12004 
iKitzbichler fc Whitd l2007l: ICirasuolo et all I2OI0I) . This 
study (and that of Henrigues et al! 2010! ) suggest otherwise 
- current models seem fully capable of reproducing the data 
given realistic assessments of population synthesis uncer- 
tainties and of the effects of observational error s. A much 
more serious problem, as already pointed out by IGuo et ah! 
l|201l! ), is that the models grow somewhat lower mass galax- 
ies too early. Objects with Af* ^ 10^° Mq are already present 
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with a large fraction of their z = Q abundance at redshifts 
of 2 or 3, whereas the observations indicate a drop in abun- 
dance by about an order of magnitude. Cosmic down-sizing 
thus appears much stronger in the real Universe than in the 
models. Reconciling theory and observation in the context 
of the ACDM cosmology will require star formation efficien- 
cies to scale with mass and redshift in a very different way 
than current models (and simulations) assume. 
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